Dynamical Correlations of the Spin-| Heisenberg XXZ Chain in a Staggered Field 
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We consider the easy-plane anisotropic spin-i Heisenberg chain in combined uniform longitudinal 
and transverse staggered magnetic fields. The low-energy limit of his model is described by the sine- 
Gordon quantum field theory. Using methods of integrable quantum field theory we determine the 
various components of the dynamical structure factor. To do so, we derive explicit expressions for all 
matrix elements of the low-energy projections of the spin operators involving at most two particles. 
We discuss applications of our results to experiments on one- dimensional quantum magnets. 
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I. INTRODUCTION 



The field-induced gap problem in anisotropic quasi 
one dimensional spin-i Heisenberg antiferromagncts has 
attracted much experimental ^"-'^^ and theoretical^' 
attention in recent years. Two scenarios have been 
studied in particular. For isotropic exchange in- 
teraction a gap can be induced by the applica- 
tion of a uniform magnetic field in presence of 
a staggered g-tensor and/or a Dzyaloshinskii-Moriya 
interaction^^. This is the case for materials such as Cop- 
per Benzoatei-6, CDC [CuCla • 2((CD3)2SO)]^ Copper- 
Pyrimidine*"^^ and Yb4As3^'^'^^. Theoretical studies 
have analyzed the excitation spectrum^^"^®, the dynam- 
ical structure factor^ '''^^j the specific heat^°, the mag- 
netic susceptibility^'^^ and the electron-spin resonance 
lineshape^^. In the materials mentioned above appli- 
cation of a uniform magnetic field H induces a stag- 
gered field perpendicular to H. It is the induced stag- 
gered field that leads to a spectral gap. The staggered 
field is generated both by a staggered g-tensor^^'"^" and 
a Dzyaloshinskii-Moriya (DM) interaction. The simplest 
Hamiltonian describing such field-induced gap systems is 
given by^^ 



-lysi (1) 



where h ~ jH . The constant 7 is given in terms of the 
staggered g-tensor^^'^^ and the DM interaction. In crit- 
ical systems with exchange anisotropy such as the spin- 
1/2 Heisenberg XXZ chain a second mechanism for in- 
ducing a gap by application of a uniform magnetic field 
exists. While application of a field perpendicular to the 
easy plane leaves the system critical, applying a field 
in the easy plane leads to the formation of a spectral 
gap 

The purpose of the present work is to extend the the- 
oretical analysis of the staggered field mechanism for 
generating a spectral gap to the case of the anisotropic 



Heisenberg chain, 

n 



(2) 



In what follows we will consider the region —1 < S < 1 
which corresponds to an "XY"-like exchange anisotropy. 
It is important for our analysis that the staggered field is 
transverse to the anisotropy whereas the magnetic field 
is along the anisotropy axis. Only in this case does the 
low-energy limit map onto an integrable model, the sine- 
Gordon quantum field theory. 

The outline of this paper is as follows: in section II 
we construct the continuum limit of the model (2). In 
section HI we derive a spectral representation of the dy- 
namical structure factor at low energies. In section IV 
we present the calculations for retarded two-point corre- 
lation functions. In section V we present our results for 
the components of the dynamical strucure factors. Sec- 
tion VI summarizes our results. The technical aspects 
of our analysis are summarized in several appendices: in 
Appendix A we discuss how the parameters of the low- 
energy field theory can be determined from the Bethe 
ansatz solution of the Heisenberg chain in a magnetic 
field. Appendices B and C present results for the form 
factors of the operators entering the calculation of the 
dynamical structure factor. 



II. CONTINUUM LIMIT 

In the limit \h\ <^ H, J the staggered field can be taken 
into account as a perturbation to the low-energy limit of 
the XXZ chain in a magnetic field. It is well-known that 
the low-energy limit of the spin-i Heisenberg XXZ chain 
with XY-like anisotropy |i5| < 1 is given by a free bosonic 
theory^^"^^ 



(3) 
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where the field ^{x) and its dual field 8(a;) are compact- 
ified 



where 



e(a;) EE e(x) + 



2n 



* = $ + Stt/?. 



(4) 



The commutation relation between $ and reads 

Q{x),^{x') ^Sm^nix^ x'), (5) 

where dnix) is the Heaviside step function, equal to 
for X < 0, 1 for X > and 1/2 for x = 0. The parameters 
V, (3 and kp (see below) in the low-energy theory can 
be calculated directly from the Bcthe ansatz solution on 
the XXZ chain'^^. How this is done is briefly reviewed in 
Appendix A. The results as well as the other parameters 
used are listed in Table I for the anisotropic parameter 
5 = 0.3. In the continuum limit the lattice spin operators 
have the following expansions 



3 

a=l 



's:{x) 



(6) 



where x = ja^ and oq is the lattice spacing. The 
wavenumbers are 



Ql = -Ql = Q = 2A:j. , 

ao 



TT 

ao 



'Ql = 2kt 



where the Fermi momentum is given by 



(7) 

(8) 
(9) 



(10) 



Here (Sj) is the magnetization per site. The continuum 
fields 5" are given in terms of the canonical boson $ and 
its dual field O as 



5f(x) 
S^ix) 

Sfix) 
5|(x) 
5|(x) 

SHx) 
S^{x) 



^A{H)[ol{x)+Olf,{x) 
c{H) cos (/3e(x)), 



^^A{H)[ol{x) ~ Olp{x) 
c{H) sin (/3e(a;)), 
Sfix)'^ 



(5|(x))^ = -la(i/)Oi(x), 



(11) 
(12) 

(13) 

(14) 
(15) 
(16) 

(17) 
(18) 



We are using normalizations such that 
(Oi(r,x)Ozi(0,0)) = 



(19) 



VT 


+ ix 




al 


VT 


— ix 




y2^2 _^ rj.2 



2a" + ^ 



(20) 

The coefficients a{H), c{H) and A{H) have been deter- 
mined numerically in Ref. [38]. The staggered magnetic 
field perturbation can be bosonized using (11) - (18), 
which leads to a sine-Gordon model 



dx 



167r 

^Ji{h,H) cos(/3e(a;)) 



{d^^{x)f + {dMx)f 



(21) 



where ^(/i, H) = hc{H). We note that as we have chosen 
to bosonize in a finite magnetic field, the cutoff of the 
theory is H rather than J. However, it is straightforward 
to recover the zero field limit (where one bosonizes at 
H ~ and the cutoff is J) in the expressions for the 
structure factor we give below. 



A. Elementary Excitations 

The sine-Gordon model is integrable and its spectrum 
and scattering matrix are known exactly^^"**^. In the 
relevant range of the parameter /3 (0 < /3 < 1) the spec- 
trum of elementary excitations consists of a soliton-anti- 
soliton doublet and several soliton-anti-soliton bound 
states called "breathers". There are altogether [1/^] 
breathers, where [x] denotes the integer part of x and 
^ = Y3gT- In order to distinguish the various single- 
particle states we introduce labels s and s for solitons and 
anti-solitons respectively and 6i, . . . , for breathers. 
Energy and momentum carried by the elementary exci- 
tations are expressed in terms of the rapidity 9 as 

wPe = Ae sinh(6') , E^^A^ cosh(6') , (22) 

where = = A, Ab^ = Ak = 2Asin(^). The 
soliton gap as a function of parameters H and h is^O-^s 



A 

7 



2v 



r(|] 



Jaoc(_ff)7r r( h 



2v 



(23) 



When (5 w 1 and the magnetization is small the leading 
irrelevant perturbation to the Gaussian model needs to 
be taken into account, leading to^^ 



A 

7 



B 



-2/3^ 



, (24) 



where B = 0.422169. 
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B. Scattering States 

It is useful to introduce creation and annihilation op- 
erators Al{6) and A^{6) for the elementary excitations. 
Here e = s, s, 61, . . . , 6[i/^]. The creation/annihilation op- 
erators fulfil the so-called Faddeev-Zamolodchikov (FZ) 
algebra 

Aai0l)Ab{02) = 5^;^>i-02)A'(02)Aa'(0l); 

AUeMlie^) = s:,'''{e,-e2)Al{e2)Ai,{e,); 
Ai{eMb{02) = sin9,-92)A,,ieMli0i) 

+ 2TTdabS{e,-e2). (25) 

Here S{9) is the scattering matrix of the sine- 
Gordon model^^"'*^. Multi-particle scattering states of 
(anti)solitons and breathers are given in terms of the FZ 
creation operators as 

\{e„,e„})^Alie„)...Alie^)\0). (26) 

Energy and momentum of these states are 

n n 

1=1 1=1 

The resolution of the identity in the normalization im- 
plied by (25) reads 

^=f:E/%#iK'M)(K,Mi. (28) 

C. Discrete Symmetries 



TABLE I: Amplitudes A, a and c, the dimensionless spin 
velocity v/ Jao, the coupling /3, and the field H as functions 
of the magnetization m for the anisotropic parameter S = 0.3. 
The amplitudes are determined in Ref. [38]. 



m 


A 


a 


c 


v/ Jao 




H/J 


0.02 


0.3044 


0.3953 


0.5275 


1.1804 


0.386192 


0.09093 


0.04 


0.3065 


0.3913 


0.5268 


1.17114 


0.385821 


0.18186 


0.06 


0.3096 


0.3867 


0.5256 


1.15828 


0.385332 


0.2598 


0.08 


0.3130 


0.3817 


0.5240 


1.13738 


0.384573 


0.35073 


0.10 


0.3173 


0.3769 


0.5219 


1.11423 


0.383768 


0.42867 


0.12 


0.3226 


0.3713 


0.5194 


1.08072 


0.38265 


0.5196 


0.14 


0.3284 


0.3661 


0.5164 


1.04600 


0.381535 


0.59754 


0.16 


0.3354 


0.3610 


0.5129 


1.01244 


0.380489 


0.66249 


0.18 


0.3433 


0.3559 


0.5088 


0.966005 


0.379084 


0.74043 


0.20 


0.3527 


0.3508 


0.5041 


0.921509 


0.377775 


0.80538 


0.22 


0.3642 


0.3460 


0.4988 


0.870927 


0.376322 


0.87033 


0.24 


0.3773 


0.3415 


0.4929 


0.813165 


0.374702 


0.93528 


0.26 


0.3923 


0.3371 


0.4861 


0.760734 


0.37326 


0.98724 


0.28 


0.4102 


0.3329 


0.4785 


0.701482 


0.371658 


1.0392 


0.30 


0.4321 


0.3286 


0.4699 


0.651491 


0.370326 


1.07817 


0.32 


0.4596 


0.3253 


0.4602 


0.575147 


0.368318 


1.13013 


0.34 


0.493 


0.3222 


0.4492 


0.507976 


0.366572 


1.16910 


0.36 


0.5342 


0.3193 


0.4367 


0.456492 


0.365244 


1.19508 


0.38 


0.588 


0.3166 


0.4222 


0.389204 


0.363518 


1.22431 


0.40 


0.664 


0.3141 


0.4053 


0.326360 


0.361913 


1.24704 


0.42 


0.769 


0.3131 


0.3851 


0.259866 


0.360218 


1.26652 


0.44 


0.934 


0.3125 


0.3602 


0.186561 


0.358349 


1.28276 


0.46 


1.214 


0.3127 


0.3279 


0.122936 


0.356722 


1.29251 


0.48 


1.89 


0.3142 


0.2796 


0.0448071 


0.354713 


1.299 



The Hamiltonian is invariant with respect to charge 
conjugation 



(29) 



The action of the charge conjugation operator C on phys- 
ical states follows from 



c\o) = |0), 
cAl{0)c-' ^ AUe), 



(30) 



We see that even breathers are invariant under charge 
conjugation, while odd breathers change sign. The topo- 
logical charge 



f3 



is a conserved quantity. We will use the conventions in 
which soliton/antisoliton and breathers have topological 
charge =Fl and zero respectively. 



III. DYNAMICAL STRUCTURE FACTOR 

The central object of our study is the inelastic neutron 
scattering intensity, which is proportional to**^ 



fc"fc" 



S"°' {uj,k). (31) 



Here a, a' = x,y, z, k denotes the component of k along 
the chain direction, and the dynamical structure factor 
on a chain with L sites is defined as 

oo 

5""'(c.,fc) = i^ / |^e-*-*'=('-'')(0|5r(i)5p'|0). 

(32) 

Substituting the low-energy expressions (6) into (32) we 
obtain 

3 ^ 
a, 6=1 

X {0\S:{t,x)SfiO,ym, (33) 



4 



where x = Iuq, y = I'ao and S"{x) are the leading terms 
in the low energy limits (6) of the lattice spin operators. 
Using that the expectation value is a slowly varying func- 
tion oi X — y we see that only terms with 

fc « Qa « -Qb (34) 

contribute to (33)'^^. The dynamical structure factor can 
be expressed by means of a Lehmann representation in 



terms of scattering states of solitons, anti-solitons and 
breathers. Inserting a complete set of states (28) between 
the operators in (33) and using 

(O|5„"(t,x)|{e„,0„}) = 

= e-'^(">*+'^<">^(O|5,"(O,O)|{e„,0„}), 

we arrive at 



5"" iLo,Q': + q) 



2tt 

ao 



n=l {e„} 6=1 



d9i... de„ 

?i!(27r)" 



'|5,"|{e„,0„})({e„,0„}|5? \Q)5{q-P{,,})5{uj-E^^}). (35) 





S! 


S2 


S3 


s! 




5| 


5! 


51 


<5| 


Q 


-1 





1 


-1 





1 





-1 


1 


c 




+ 

















TABLE II: Topological charge Q and eigenvalue (where ap- 
plicable) under charge conjugation C of the continuum spin 
operators. 



A 



Here q is assumed to be sufficiently small {q 
kp, ^). Due to energy-momentum conservation only a 
finite number of intermediate states contributes to the 
correlator (35). Moreover, at low energies contributions 
of intermediate states with large numbers of particles to 
the correlator (35) are generally small^^'^^. We therefore 
restrict our following analysis to one and two particle 
contributions. Many matrix elements in (35) are in fact 
zero as can be established by using charge conjugation 
symmetry and topological charge conservation. The rel- 
evant properties of the continuum spin operators 5" are 
summarized in Table II. Using these properties we fur- 
thermore conclude that at low energies the non- vanishing 
components of the dynamical structure factor are 

1. S^^, in the vicinity of the points k = ±Q; 



2. S*^^, near the point fc = ^; 

3. S^^ in the vicinity of the point fc = 0; 

4. S*^^ near k = ±2kF. 

In the following we determine these in the "two-particle 
approximation", i.e., keeping only terms with n < 2 in 
the spectral representation (35). In order to do so we 
make use of the exact form of the matrix elements enter- 
ing the Lehmann representation, which follow from the 
form-factor bootstrap approach^°'^^. 

We note that as a consequence of charge conjugation 
symmetry the components of the structure factor in the 
vicinities of fc = Qa and k = —Qa are the same. 

IV. CALCULATION OF CORRELATION 
FUNCTIONS: KINEMATICS 

The formalism we employ to calculate the dynamical 
structure factor can be used quite generally to determine 
(real and imaginary parts of) two-point correlation func- 
tions. The retarded two-point function of two bosonic 
operators A and B has a spectral representations of the 
form 



n=l{e„} 



d9i... de„ 



(0|A|{e„, 6'„})({ 



5{vq - P{„}) 



n!(27r)" 

{0\B\{ })({£„, 0„}|A|O) 
I 



5{vq + P{^^) 



^ + £^{n} + if] 



(36) 



Here 77 is a positive infinitesimal, |{e,i, On}) are n-particle (22) respectively. The leading contribution to the spec- 
scattering states of solitons, antisolitons and breathers tral sum in (36) is due to intermediate states with one 
(26) with energies and momenta arc given by (27) and and two particles. Using momentum conservation it is 



5 



possible to simplify the expressions for these contribu- butions to G^^ 
tions as we discuss next. 



A. One-particle kinematics 

Resolving the momentum conservation delta function 
leads to the following result for the one-particle contri- 



J 



E 



aoeaiq) 



U! — Aa cosh + 11] 



uj + /S.a cosh6' + iri 



^ - £a{q) + irj 



LO + £a{q) + irj 
I 



(37) 



where a runs over all single-particle labels (i.e. soliton, 
antisoliton and breathers) and 



ea{q) - ^A2+w2q2 
6'q = arcsinhf— ^ 



B. Two-particle kinematics 



As two-particle form factors of scalar operators depend 

(38) 

only on the rapidity difference, it is useful to change vari- 
^3g^ ables to 9± — {6i ± 6*2)72. Resolving the momentum 
conservation delta function then gives 



V ^ r de 1(102 



{0\A\92,9i)a2ai aia2(6'li ^2|-B| 



^{""q - Z]j=i sinhflj ) 
uj - Y.]=i cosh6'j + irj 
5{vq + ^a, &\-ah6j) 



V y-^ r de_ 

a,b 



\B\92,0l)a2ai aia2{(^li(^2\A\0) 

(OlAie^" - 9^,9f + 9_),a ab{et + 9^,9f - 9_\B\0) 



^ + Ej=l C0Sh6'j + IT], 



eab{q,9-) {lJ - eab{q,9-) + 17]) 

{0\B\9S''- 9-,9f + 9^)ba ab{9f + 9^,9^'' -9^\A\0) 
eabi.q-,9-) {uj + eab{q,9-) + iij) 



(40) 



where 



where 



£ab{q,6) 



v^q^ 



A2 + A2 + 2A,Ab cosh(20) 



9^'' = In 



vq + £ab{q,9^) 



Aaexp(6i_) -t- Abexp(-6i_) 
The imaginary part of G2^ {oj, q) can be simplified using 
1. 1 



-Im - 

TT Sabiq, 9^) [uJ - eab[q, 9-) + 11]] 

_ 5{9^-^^) + 5{9^ + ^^) 



v/s2 - (A, + A,) - (A, - A,)- 



F, (41) 



2 2 

V q , 



9ab{s) ~ arccosh 



A?,-A2i 



2A„A6 



(42) 
(43) 



Carrying out the 0_ integral using the delta functions we 
obtain 



6 



-IinG^p^(a.>0,g) 



2nai 



■EE 

a,b <y=± 



{Q\A\eL{^. q)ba ab{ej{u^. q),OLi^, q)\B\o) 



(A, + A,) " (Aa - A,)2 



(44) 



where i?_f/(x) is the Heaviside function and 

1 



9^i^{uj,q) — arcsinh 



2A,s2 



vqis' + Ai- At) 



± LU^is^ - (A, - A,)2)(s2 - (A, + Afc)2) 



(45) 



and five breathers with gaps 

Ai = 0.54098A , A2 = 1.04162A , A3 = 1.46461A, 
A4 = 1.77841A , A5 = 1.95962A. (47) 

In order to broaden delta functions appearing in one par- 
ticle contributions, wc introduce a small imaginary part 
in w, equal to ry — O.OIA. 



The two terms in (44) arise from the two delta functions 
in (41). Using the results summarized in this section 
we can determine the one and two particle contributions 
to both real and imaginary parts of two point functions. 
The two particle contributions to the real part involve one 
(principal part) integration, which is readily performed 
numerically. In order to determine the dynamical struc- 
ture factor we only require the imaginary part of several 
two point functions. 



V. RESULTS FOR THE DYNAMICAL 
STRUCTURE FACTOR 

Below we present results for the dynamical structure 
factor 5'"^(w, + q) (35) in the regime -1 < 5 < 1 and 
for magnetic fields H < He = J{1 + 5). Wc note that if 
H He or 6 ^ —I the cut-off in the field theory is very 
small, which limits the utility of our approach. For the 
sake of clarity we use a particular set of parameters in all 
plots 



7 = 4= 0.01191, 
H 



S = 0.3, H = 0.2598 J. 



(46) 



These correspond to a magnetization per site of (S^) — 
0.06 (see Table I) and ( = 0.174371. The spectrum con- 
sists of soliton and antisoliton with gap A w 0.04897J 



A. ^""(tJ,/;) 

In the continuum limit S^^{uj,k) is non- vanishing in 
the vicinity of the points k = ibQ and Wc will con- 
sider both cases in turn. As we have noted before, the 
response at fc = ±Q is indcntical as a result of charge 
conjugation symmetry, so that it is sufficient to consider 
k^-Q. 



1. Momenta k ^ -Q = IS') 

In the continuum limit S'^'^{uj,—Q + q) with q <^ Q 
is given by the two-point function of 5f with Sf (11). 
This is because vQ is a large energy scale proportional to 
the cutoff in the theory. Using Table II wc find that the 
following intermediate states with at most two particles 
contribute 

1. Single-soliton states. 

2. Two particle states containing one soliton and one 
breather. 

The corresponding matrix elements are calculated in Ap- 
pendix B. Using the results of section IV to carry out the 
rapidity integrals we arrive at the following expression for 
S^^{uj, —Q + q) within the two-particle approximation 



S'^^'ico > 0,-Q + q) 



fe=i 



V(,s2-(A-Afc)2)(,s2-(A + Afc)2) 

A7.fMs..(s))e-''-^(-''' " 



cr=± 



(48) 



r 



Here s is the Mandelstam variable (42), the overall nor- where Zi{(3) is given by equation (B2), the minimal 
malization is form factors i^s™™(^) by equation (B6), the pole func- 

1/2 



A = A{H)(zi{(3) 



(49) 
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tions Kff^^ (0) by equation (B22) for k even and (B23) for where 
k odd, the normahzation factor N^^^^ by equation (B20) 
and the functions 6sb^ (s) and O'^^^^ (w, q) arc presented in 
equations (43) and (45), respectively. 



0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 



y/A 



FIG. 1: One and two-particle contributions to S^^{lj, —Q + 
A/v) as a function of uj for S = 0.3 and H = 0.2598 J. The 
delta-function peak (pink) has been broadened to make it 
visible. 



2(/3 



4/3 

'2"^ while 



)^ > 1. For large oj this in- 
is goes to zero in a power-law 



creases as oj 

fashion for lo vq. In presence of a staggered field, 
the dynamical structure factor (48) has divergence for 
u; ^ ^(A + Afc)2+i;2g2 (fc = 1, 2, . . . , [1/^]), while the 
large frequency behavior is the same as without the stag- 
gered field. 



2. Vicinity of antiferromagnetic wave number: k ^ ir/ao 



In the continuum limit 



q) with qaQ <C tt is 



given by the two-point function of the charge neutral op- 
erator iSf (13). Using Table II and (30) we find that the 
following intermediate states with at most two particles 
contribute to the two-point function of 5f 



We note that S^^{uj, —Q + q) vanishes when q ^ 0. In 
Figure 1 we therefore plot S^^{uj, —Q+A/v) as a function 
of w. In order to broaden delta-function contributions we 
introduce a small imaginary part in lo. Two features are 
clearly visible: there is a coherent peak corresponding 
to the contribution of single-soliton excitations at energy 
at AV2. At higher energies breather-soliton continua ap- 
pear. Their contributions grow with increasing lu because 
5f is an irrelevant operator. It is instructive to compare 
our result to the gapless spin-1/2 Heisenberg XXZ chain, 
see e.g. Rcf. [52]. There one has 

5-(^ > 0, -Q + g) (X f '^f,^' , (50) 
I 



1. Single breather states even under charge conjuga- 
tion, i.e. BUom. 



2. Two particle states containing one soliton and one 
antisoliton. 

3. Two particle states containing two even or two odd 
breathers. 



Using the results of section IV, we obtain the following 
expression in the two-particle approximation 



5'^^ ( a; > 0, — + g 
ao 



TTao 



fe=l 



,(s))rz?ff(s-2A) 



[i/«] 

V St 



rcvcn I p, 
/ . ''k+k'l^bkb 
k,k' = l 



(^fc.6„(.s))| 



4A2 

-AfeO 



V(s2 - (A,. - Afe02)(s2 - (Afe + AfcO^) 



(51) 



Here the single-breather form factors Fj^^ are given 
by equation (C14), the soliton-antisoliton form factor 
F^°^'"^^\0) by (CIO) and the breather-breather form fac- 
tors fII^^^X^) by (C24) respectively. The function e,,.{s) 
is given by (43) and 

jovon ^1 1, if fc is even, ^^^^ 
^ 1 0, overwise. 

In Fig. 2 we plot the dynamical structure factor (51) 
as a function of frequency. We note that because iSf 



is a scalar operator S'^^(a', ^ + ?) depends only on the 
Mandelstam variable s (42) rather than on uj and q sep- 
arately. The first peak in ^^'^(ti;,^) is due to the 62 
singlc-brcather excitation (blue line). At = A4 there is 
a second single-breather contribution, due to 64. Above 
LO = 2Ai a strong 6161 two-breather continuum occurs 
(pink line). Around lo = 2 A contributions from soliton- 
antisoliton and 6163 and 6262 two-breather continua are 
visible. We note that the thresholds of 6163, ss and 6262 
continua all occur around 2A is a peculiarity of the pa- 
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a2 I 




FIG. 2: One and two-particle contributions to S^^{ijj, ^) as 
a function of oi for 5 = 0.3 and H — 0.2598J. Delta-function 
peaks (blue) have been broadened to make them visible. 

rameters we have chosen in the plots. 

B. S'f«(ix),fc) 

Next we turn to the j/y-component of the dynamical 
structure factor. In the continuum limit S'^^(tJ, k) is non- 



vanishing in the vicinity of the points k = ±Q and 
We will consider both cases in turn. 



1. Momenta k -Q = ~^{Sj) 

In the continuum limit S"^^ {to, + q) with q <^ Q 
is given by the two-point function of iS| with Sf (14). 
Using Table II we find that the following intermediate 
states with at most two particles contribute to the two- 
point function 

1. Single-sohton states. 

2. Two particle states containing one soliton and one 
breather. 

The corresponding matrix elements are calculated in 
Appendix B. Carrying out the rapidity integrals, see 
section IV, we arrive at the following expression for 
S^y{ijj, — Q + q) within the two-particle approximation 



syy{uj > 0,- 



q) 



2vA'^ v'^q^ + 



Sis' - 



^2 If™" 



47rao 



''"^ V(s2_(A-Afc)2)(s2-(A + Afe)2) 



X E [i^f,J^^sbJ.))e^-'=(--^)-iff,J-a0.,,(s))e''^-.(-'' 

I 



cr=± 



(53) 



Here the overall normalization A is given by equation 
(49), the minimal form factor F^l'^ie) by (B6), the pole 
function i^rfbjei) by (B22) for k even and (B23) for k 
odd, the functions Ogb^ (s) and 0^^^ {to, q) by (43) and (45), 
respectively. 



12 




0.004 
0.003 






10 
8 




0.002 
0.001 






6 






1 


2 3 4 


4 










2 












1 


2 


3 


4 5 



co/A 



FIG. 3: One and two-particle contributions to S^''{uj, —Q + 
A/v) as a function of a; for 5 = 0.3 and H = 0.2598 J. The 
delta-function peak (pink) has been broadened to make it 
visible. Inset: the soliton-breather two-particle contributions. 



make them visible. We see that there is a coherent peak 
corresponding to the contribution of single-soliton exci- 
tations at energy at At higher energies breather- 
soliton continua appear. Their contributions grow with 
increasing lu because iS^ is an irrelevant operator. 



2. Vicinity of antiferromagnetic wave number: k ~ Tv/ao 

In the continuum limit Sy^luj, + q) with qaQ <C tt is 
given by the two-point function of the charge neutral op- 
erator iSf (15). Using Table II and (30) we find that the 
following intermediate states with at most two particles 
contribute to the two-point function of 5| 

1. Single breather states odd under charge conjuga- 
tion, i.e. BL+i(e)|0). 

2. Two particle states containing one soliton and one 
antisoliton. 



We plot S'^y{uj, — Q-f A/ti) as a function of w in Figure 
3. Delta- function contributions have been broadened to 



3. Two particle states containing one even and one 
odd breather. 
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Using the results of section IV, we obtain the following 
expression in the two-particle approximation 



[uj>o, — + q 

ao 



vc^H) 



k=l 



Ft 



sVs^ - 4A2 



fc,fc'=i 



- (A, - Ak'r){s^ - (Afc + AfcO') 



(54) 



Here the single-breather form factors are given by 



(C14), the soliton antisoliton form factor F, 



sin(/3e) 



(0) by 



(Cll), the two-breather form factors -F;,b,,(^) by (C24), 
the function Oab{s) by (43) and 



4 



dd _ I 1, if fc is odd, 
I 0, overwise. 



(55) 
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FIG. 4: One and two-particle contributions to S^''{uj,-^) 
as a function of u for 5 = 0.3 and H — 0.2598J. Delta- 
function peaks (blue) have been broadened to make them 
visible. Insert; the soliton-antisoliton and breather-breather 
two-particle contributions. 

We plot Syy{uj, ^) as a function of uj in Fig. 4. We 
see that it is dominated by the contribution of the first 
breather bi (the corresponding delta function has been 
broadened). The contributions from 63 and 65 single- 



breather states are small in comparison. Similarly, the 
two-particle ^'1627 ss and 6164 continua shown in the inset 
of Fig. 4 are negligible. 



C. Longitudinal structure factor S^^{uj,k) 

We now consider the zz-component of dynamical struc- 
ture factor. In the continual limit S^^lujjk) is non- 
vanishing in the vicinity of the points fc = and ±2kp. 
We will consider both cases in turn. 

1. Vicinity of ferromagnetic wave number: fc ~ 

In the continuum limit S^^{u>^ q) with qao <C tt is given 
by the two-point function of the charge neutral opera- 
tor Sf (17). Using Table II and (30) we find that the 
following intermediate states with at most two particles 
contribute to the two-point function of Sf 

1. Single breather states odd under charge conjuga- 
tion, i.e. sL+iWlO). 

2. Two particle states containing one soliton and one 
antisoliton. 

3. Two particle states containing one even and one 
odd breather. 

Using the results of section IV, we obtain the following 
expression in the two-particle approximation 



S^\io>0,q) 



k=l 



2 a,b^u;^ \F^MM)\^h{s - 2 A) 
^[s -Ak) + 



aob'^Lo'^ 



E '5^+1-' In!.,, (^....(-)) I 
I — 



sVs"^ - 4A2 
dnis-Ak-Ak' 



k,k' = l 



V(s2 - (Afc - A,,y){s^ - (Afc + Afc,)2) ' 



(56) 



where the single-breather form factor F^^ is given by by equation (C12), the breather-breather form factor 
equation (C15), the soliton-antisoliton form factor F^{9) 
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^bkbyi^) equation (C26), 9ab{s) is given by equation 
(43), s is the Mandelstam variable (42), 5'^'^'^ is given in 
(55) and the overall normalization is 




FIG. 5: One and two particle contribution to S~~{uj,G) (56) 
as a function of oj for 5 = 0.3 and H = 0.2598J. Delta- 
function peaks (blue) have been broadened to make them 
visible. Insert: the soliton-antisoliton and breather-breather 
two-particle contributions. 

The dynamical structure factor (56) is shown in Figure 
5. Note that since SI is a scalar operator, S^^{uj,q) de- 



pends on the Mandelstam variable s (42) rather than on 
uj and q separately. In order to broaden the delta function 
contributions we introduce a small imaginary part in lu. 
The dominant peak in S^^{uj,q) is due to a bi breather 
contribution. The contributions due to 63 and 65 breather 
states are much smaller. The soliton-antisoliton and 
breather-breather contributions to 5*^^(0;, q) are barely 
visible in the figure. 



2. Momenta fc ~ — 2fcF 

In the continuum limit S^^'{lo, —2kp + q) with qao <^ tt 
is given by the two-point function of 5| with 5| (18). 
Using Table II and (30) we find that the following inter- 
mediate states with at most two particles contribute to 
the two-point function of 

1. Single soliton state. 

2. Two particle states containing one soliton and one 
breather. 

Using the results of section IV, we obtain the following 
expression in the two-particle approximation 



S'-'^ico > 0, -2kF + q)~ S{s^ 

ao 



~2 2 



r 



A,) 



v/(s2-(A-Afc)2)(s2-(A + Afc)2) 



(58) 



Here the minimal form factor F^l™ 

sou 



[0) is given by (B6), 
the pole function K^bki^) ^y (B22) for k even and (B23) 
for k odd, the function Bsb^is) by (43), the overall nor- 
malization is 



a(H] 



Z{0) 



(59) 



The dynamical structure factor (58) is shown in Figure 
6. Here we chose q = 0. The strong low-energy peak 
in S^^ {lli, —2kp) is due to a one-soliton state. Soliton- 
brcather continua appear at higher energies. 



VI. SUMMARY AND CONCLUSIONS 

In this work we have determined the low energy dy- 
namical spin response of the anisotropic spin-1/2 Heisen- 
berg XXZ chain in the presence of both uniform and 
stagcred magnetic fields. The uniform field was taken 
to be along the anisotropy axis and the staggered field 



perpendicular to it. The qualitative features of the model 
such as a field induced gap and the formation of bound 
states are similar to the case of isotropic exchange, which 
has been previously studied in detail^'^^"^*'^''"^^'^'''^^. 
The main effect of a strong exchange anisotropy is to 
generate further bound states and increase the bind- 
ing energy. We have analyzed these effects on the dy- 
namic response and determined for the first time all 
two-particle contributions, in particular those contain- 
ing one soliton and one breather. The results obtained 
here can be used to study a quasi one dimensional array 
of anisotropic Heisenberg chains in a uniform magnetic 
field by combining a mean-field approach with an RPA- 
like approximation^'^^^'^ . This is of interest in view of neu- 
tron scattering experiments on the quasi-lD anisotropic 
Heisenberg magnet Cs2CoCl4^2. 
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FIG. 6: One and two particle contribution to S^^ {ui, —2kF) 
(58) as a function of lo for S = 0.3 and H = 0.2598J. Delta- 
function peak (blue) has been broadened to make it visible. 
Insert; the soliton-breather two particle contributions. 
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APPENDIX A: SPIN VELOCITY, FERMI 
MOMENTUM AND COMPACTIFICATION 
RADIUS 



We note that spinons can only be excited in pairs. The 
magnetization per site in the ground state is given in 
terms of the ground state root density p{X) as 



dX p{X) 



The Fermi momentum is equal to 
2n 



kp = p{A) = 



dX p{X) 



ao 



-is;) 



(A4) 



(A5) 



where we have used that p(— A) = p{X). The spin veloc- 
ity is equal to the derivative of the spinon energy with 
respect to the momentum at the Fermi points 



deiX) 



dp{X) 



de{X)/dX 



ao 



x=A 27rp(A) 
Finally, the dressed charge is related to /? by 

,3^ ' 



^/^Z{A) 



(A6) 



(A7) 



In order to determine v and (3 wc solve (Al) numerically, 
which is easily done to very high precision as the equa- 
tions arc linear. The results are shown in Fig. 7, 8 and 9. 



In this appendix we summarize how to determine the 
parameters of the Gaussian model (3), (4) that describes 
the continuum limit of the Heisenberg XXZ chain in a 
magnetic field from the Bethe ansatz solution^^. The 
velocity, Fermi momentum and compatification radius 
are expressed in terms of the solutions of the follow- 
ing set of linear integral equations for the dressed energy 
e(A), dressed momentum p(A), dressed density p{X) and 
dressed charge Z{X) 



e{X) 
Z{X) 



^ KiX - p) e{p) = H- 



J sm 7 



cosh2A — cos 7 



— I dp p{p) 
ao Jo 



dp 



2 sin 7 



_^^^K{X-p)p{p) 2.[^„^i^2A-cos7] 

A 



I ^K{X~p)Z{p)^l 
J~A 27r 



(Al) 



Here the exchange anisotropy is parametrized as S 
cos(7) and the integral kernel is given by 



A'(A) = -2sin27/(cosh2A - cos27). 



(A2) 



The integration boundary A is fixed by the condition 

ei±A) = . (A3) 

The physical meaning of the various quantities is as fol- 
lows: e(A) and p(A) are the energy and momentum of an 
elementary "spinon" excitation carrying spin = ±i. 



8=0 

8=0.5 

8=1 




FIG. 7: Spin velocity as a function of magnetic field for dif- 
ferent values of S 

For zero magnetic field we have kp = 7r/2ao, /3'^ = 
arccos(— 5) and 

JflQ sin47r/3^ 



2 1 - 4/?2 



(A8) 



APPENDIX B: FORM FACTORS OF THE 
EXPONENTIAL FIELD O^^ 

The exponential field (19) has topological charge 
4-1. Hence it has non-vanising matrix elements between 
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nil ^, 



5=0 
8=0.5 
5=1.0 




the ground state |0) and one soliton states, two-particle 
soliton-breather s6„ states, etc. In this appendix we de- 
termine all one and two-particle form factors. Our results 
hold generally for the sine-Gordon model in the attrac- 
tive regime. 



FIG. 8: "Fermi momentum" kp as a function of magnetic 
field for different values of 6 



1. One soliton form factors 



0.5 



0.45 



0.4 



0.5 



1 

H/J 



5=0 

8=0.5 

8=1.0 



1.5 



The one-soliton form factor is^^ 



(Bl) 



FIG. 9: Parameter /? as a function of magnetic field for dif- 
ferent values of 5 where 



J 



1/4 



2Cf 



X cxp ■ 



16 



-1/8 



V^Aaor((3 + 0/2) 



-| 2a^ + l/8/3^ 



dt 



cosh(4Cat//3)e-(i+«)* - 1 



t V4sinh(t^) sinh(t(^ + 1)) cosh(t) 1 siiih(t^) 



-2t 



(B2) 



2. S-matrices, their analytical properties and 
minimal form factors 



The soliton-breather S-matrix is given by' 



57 
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S,b„iO) = (-l)"exp<^2 



d< cosh(^t) sinh«t) . ^ f 2dt 
t cosh(t) sinh(^t) 



1 sinh ( 9 



n-g(n+l-23) 



r 



(B3) 



The corresponding soliton-breather 3-particle coupling is where TZsb„ {S) is given by 

1/2 



9sb„ 



2o„.(=|i)n».»(^ 



1=1 



Note that by crossing symmetry we must have 

s _ b„ 

9sbn — 9SS ■ 



(B4) 



(B5) 



The minimal soliton-breather form factor can be ob- 
tained combining equations (2.23), (4.19), and (4.20) of 
Ref. 51, which can be summarized as 



Rsb2„ (^) — 1 1 

^..&2„+i(^) = ismh{e/2). (B7) 



3. Soliton-breather form factors 



We will calculate the full form factor using the residue 
condition'"'^'^* 



^(e) exp <; / dtf{t) sinh [ — ]}, 



i^"""(^) =exp<^ / dtfit) 



1 - cosh[t(l + A)] 
2 sinh(i) 



The soliton-breather S-matrix (B3) then gives rise to the 
minimal form factor 



FZ\e) = 7^.,„(0)cxp 



dt cosh(^t) sinh(n^i) 



t cosh(t) sinh(^t) 
^ 

1 - cosh[2t(l + £)] \ 
sinh(2t) I ' 



WssF}CAOs,0b) 



Ress=oFt,{6 



where u^^ = 7r(l — n^), and 



9-sl = (-1) 5J 



n b„ 



2 cot 



2 



n 

1=1 



coV 



~2~ 



(B8) 



1/2 



(B9) 



The three particle form factor involving two solitons and 
(B6) one antisoliton is^^ 



71, (72, (73 



) = 



iC2^yZl{a) V 



d0o 



4Ci 



/ ^e-(l?+t)»«W7(0io)VF(02o)VF(03o)-e^ / ^e~'-^^^^>'"'W{eio)W{92o)W{9o3) 



e3]ret(«i+»^+«3)+f G(0i2)G(0i3)G(023) 



where 9jk = 9j — ' 



2tt 



G{9) = iCi sinh ( - ) exp <; - 



d9 



2tt 



C- 



w{e) 



cosh(6i) 



exp < — 2 



dtsinh((l~g)t) sinh"(t(l - f )) 
t sinh(2t) cosh(t) sinh(^<) 



dt sinh(i(C - 1)) sinh^(t(l - i9/Ti)) ] 
T sinh(2i) sinh(t^) I 



, (BIO) 

(Bll) 
(B12) 



r 



Ci = exp 



C2 = exp 



dtsinh^(</2)sinh(<(C- 1)) 
T sinh(2t) sinh(t^) cosh(t) 

'dtsinh^(t/2)sinh(t(^- 1)) 
1 sinh(2t) sinh(t^) 



The integration contours (7+ and C_ are constructed 
(B13) as follows. The contour C+ runs from —00 to 00 in the 



(B14) 
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complex 0Q plane, passing above the poles at 9p + i7r/2, 
p = 1,2,3. Similarly, the contour C_ runs above the 
points 9p + iTr/2, p ^ 1,2, and then below 63 ~ in/ 2 (see 
figure 10). 




FIG. 10: Integration contours C+ (panel (a)) and C- (panel 
(b)). Both the contours are located in the strip — 7r/2 — < 
•Q60 < 7r/2 + 0. 

Let us consider the form factor (BIO) for rapidities 

01 = 9,, 62 = Ob - and 9^= 9^ + 6+ The 
function W{9) has poles in the strip |Im 6'| < -I of the 
complex 0-plane at the points 9 = ivk, where 



0<k < 



As a result the function W {9io)W {92o)W {9o3) has poles 
at 6*0 = 6*1 - iTT{k^ - i), 9o = 9b + ^(n - 2k), and 
6*0 = 9b + S-'^{n-2k), 0<k < f + By construction 
the contour C_ runs between n + 1 pairs of poles at 9q = 
9b+'^in-2k), and 9o = 9b + S+'^{n-2k), < k < n, 
with only an infinitesimal separation S between them (see 
Fig. 11). As a result the integral over 9q exhibits a 
simple pole for S 0. In order to extract the residue 
of this pole, we deform the "singular" contour C_ into a 
"regular" contour C plus closed contours including one 
pole from each pair. C is chosen such that the integral 
over it is finite in the limit (5—^0. The contours C and 
C- are shown for n = 1 in Fig. 11. 
For general n we then find that 



^e-(t+|)»oM/(0io)M^(^2o)M^(^o3) 



W{9 



lb 



i7r^(n — 2k) ' 



X Wrcs{ivk)W{-6 + ivn-k) + regular part, (B15) 



where Wros(wfc) denotes the residue of the function W{9) 
taken at the point 9 = ivk- We note that W{—S + iVn-k) 
has a simple pole at 6 = 0. 

The regular part of the integral does not contribute to 
the residue in the right hand side of equation (B8) and 
in the following will be ignored. Similarly, the integral 
over C+ has a finite limit 6^0. 
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FIG. 11: Integration contour C- for ^ = 0.174371 and n = l. 
The contour C- is transformed into the contour C and the 
closed contours around the poles. 

In order to proceed we need to analytically continue the 
function W{9), which can be done using the rclation^^ 



W{9 ± iirC) 



sin [1(0 Tf)] 



-W{9). (B16) 



sin [^{9±iTii±^-)\ 
Using (B16) the residues of W{9) arc readily calculated 



fc-i 



VC2Sm(-2^) \ ^ 



A: > 1, 



ICS {iV(^ 



/Co 



(B17) 



The soliton-breather form factors can now be deter- 
mined from the residue condition (B8) 



X G 9sb + 



^b) 

iui 



iC2\JZi{a) 



fc=0 



2 

-k. 



G 



4Ci 



^sb 



eW+n^Gi-iuf) 



cs 



X W[ 



7s& 



.<(n- 2fc) 



(B18) 



The soliton-breather form factor can be expressed in 
terms of the minimal form factor (B6) as 



FltiOs^Ob) = Ntb^^z;{p)e'w+^K%^{9sb)x 

xFt:{9sb). (B19) 
Here TV^j^ is a normalization constant given by 

(B20) 
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where 



,sbi. 



= 77(^e + l)- 



(B21) 



The "pole functions" K^i^ {9) have somewhat different 
forms for even and odd n respectively and are given by 



m — 1 



cosh 



+ eosh(f - 



1 cosh 



k=l 



1 1 

n ^^^h(f + Hiljli) cosh( 

m 

X ^ Wros(«Wm-fc)VFros("'m+A;+l 



2 2 

Ta(2fc+l)i; 



/a 



fc=0 



where 



ixa(2fc + l)i; 

«'2fc+l(6') + e 3 ^'2fc+i(-^ 



r 



(B22) 



(B23) 



2" ^ ^ 



, (B24) 



Evaluating the integrals gives 



Sbkbi{0) = tanh 



coth 



e - iu' 



bb 
l+k 



fie , 

^2 + 2 ) 



sm 



,66 



(B25) 



APPENDIX C: FORM FACTORS OF exp[iae] 



X tanh 
X J]^ < tanh^ 

X coth^ 



2 



/ — k \ I ^ — k 

coth 



,66 

^l-k+2j 



(C3) 



The operator exp[mO] carries zero topological charge. 
Hence the non- vanishing form factors with less than three where now k < I, and 
particles involve a single breather, a soliton-antisoliton 
pair or two breathers. In this appendix we construct 
these one and two particle form factors. 



,66 



(C4) 



1. S-matrices and minimal form factors 

The breather-breather S-matrix is^^ 
dt sinhm) 



The simple pole in Sb,.bi{0) at 6* ~ corresponds 
to the formation of a b^+i breather bound state. The 
corresponding residue is 



SU^^,) = 2.tan(i^iM)tan(M)cot(^) 



{CXJ 
4 J 
n 



t cosh(t) sinh(^t) 



(CI) 



k-l / ■ /■\ l+k-l 



X cosh(^t) sinh(fc^t) cosh((l - lS)t) | , fc < Z, 
^ dt sinh(f^) 



The corresponding three-particle coupling is 



^fc.fc.W = -cxp<^2 



t cosh(t) sinh(f t) 



(C2) 



bk+l 
Sbkbi 



2 cot (M) tan (M Van ^(^^^^'^^ 



cosh(^i) sinh((2fc^ - l)t) + sinh((l - $)t) 



X n 



cot , tan ( (i±|M 



(C6) 
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The minimal form factor for two different breathers is are^^ 
then 



dt cosh(Ci) cosh((l - im ^cos(p0),as _ G0 G{0) /ttA 8^cosh(f) 



X cosh 

(C7) V 2^ 







sinh(fcCi) cosh(2t(l - f )) - 1 1 X cosh ( ^ ) , (CIO) 



sinh(^t) sinh(2t) j 

The minimal form factor involving two 

breathers of the same type is given by 



(o\ j J^< ..ni.'wi-f)) f;;*«)(<,) = -|S^c„t(f)Ji=!gL 

^ /9~iTT\ 



cosh(^i) sinh((2fcC - l)i) + sinh((l - ^)t) I " V 2^ 



X 



xsinhl— (CU) 



cosh(t) sinh(^t) 



(C8) 



2. Soliton-antisoliton and one-breather form factors 

The soliton-antisoliton form factors for the operators 
cos(/3e), sin(/3e), and 9, 

za© 



6 = -i hm 9ae*°^, (C9) respectively, and 



Ci /3cosh(2-2|2^) cosh(|) 



where G{9) and Ci are given by eqs. (Bll) and (B13), 



-I 2a^ 



v2r{ 



/32 



2-2/3^ - 
■) 



2-2^: 



dt 



- 20^6-2* 



sinh^(2a/3t) 



2sinh(t/32)sinh(t)cosh(t(l - 



(C13) 



The single-particle form factors and for the oper- 
ators e*''® and Q, respectively, can be obtained from the 
residue condition for the soliton-antisoliton form factor 
and FfM 

g'^-.F.l = Rese=-»- < = 7r(l - <), 



where a = (3,&. Using equations (CIO), (Cll), and 
(C12), we can write 



2cot(^^^ sin(riO(-i)" 

2cot(^)nr->ot2(^) 



—Gp cxp 



sinh(t(^ - 1)) sinh^(tO- 
T sinh(2t) cosh(t) sinh(t^) - 



z(-l)' 



_iG(*7r((2m-l)e-l)) 



/5 Sin ((^12^) 



(C14) 
(C15) 



3. Breather-breather form factors of exp(ia9) form factor Fj^^ (7) from Ref 



60 



To calculate two breather form factor F^^^ (^12), we F"" (^] ^ K"'('A TT R('y (Tlfi'l 

will start with the formulas for n-breather (n = fc + ' 2»/2 "^-^^ ^ ^ 
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Ril) 



sinh(i(7 - zTT^)) sinh(i(7 + iO) 



, (C17) 



^«(7) = E---E(-i)''^'"^'"Ka)x 

;i=o i„=o 



X n 

l<i<j <n 

2 



l + (/»-?j)- 



i sin(7r^) 



sinh(7jj) 



(C18) 



itt) sin(7r^) 



n/2 n 



n 

i=l 



e /3 



-(-1)' 



(C19) 



where 7 = (71, . . . ,7„), / = (Zi, . . .,/„), 7^^ = 7i - 7i> 



m = 1, 2, . . . , fc, and are infinitesimal parameters. 

It should be noted that the order of calculation of 
residues predicts the rules 

l^il < |'52| < ••• < < I^i2lwc22) 

|ei| « lesi < ••• « « 1^121. 

The n-particle (n = fc + /) form factor (C16) depends 

on 7i — 7j, 1 < i < J < n. Taking into account equations 

(C21), (C20), and (C22), we can write 

-ffc-j - i{i - if («, j) e Ai, 

7,;^- = -5i^j - i{j - i)TT^, if (i, j) e yl2, (C23) 



JV = - exp 2 



(itsinh(^t) sinh((l-^)t) 
1 cosh(t) sinh(2i) 



The normalization factor Ga/'2^^^ is introduce in Ref. 



59 



Form factor Fi 



bkb. 



92) can be calculated as follows. 



• First step, we calculate the residues of the form 
factor at the point 7„_i — 7n = — «7r^ iln-i ~ 

in-i - In = 7^-1 + ^1 + As a resuh, 

we obtain the form factor for n — 2 breathers bi 
and one breather 62- 

• Second step, we calculate the residue of the fc+Z — 1- 
particle form factor obtained at the point 7„_2 ^ 

7^.-1 = ^ (7n-2 = 7^1-2 - in~l = 7^i-2 + 

82 + As a result, we obtain the form factor for 
n — 3 breathers bi and one breather 63. 

• Step number Z— 1, we calculate the residue of the k+ 
2-particle form factor obtained at the point 7a:+i — 

7fc+2 = (7fc+i = 7fc+i- 2 ' 7fc+2 = 7fc+i + 
bi-\ + ■^). As a result, we obtain the form factor 
for fc breathers b\ and one breather bi. 

Finally, taking 7^_|_i = ^2, we can write 

7fe+m = ^2 + 2^ ^ ,(C20) 

J— / — ?n+l 

m = 1,2, and bj are infinitesimal parameters. 

Similar calculations performed with the variables 7^, 
j = 1, 2, . . . , /c, give 



fc-i 

_7 — fc — 



- '(^'^":+^^"^ (C21) 



where the manifolds ^1,2,3 are constructed as following, 



^1 
^2 
^3 



1 < i < j < fc, 

fc + l<i<j<'^, 

1 < i < fc, fc + 1 < j < n. 



Then we obtain the following expression for the two- 
particle form factor. 



K6,^6".6,(^i2)^^a"(^i2), (C24) 



where the minimal form factor F^^^^(d\2) is given by 
equation (C7) for k < I, and equation (C8) for k = I, 



n 



sinh (i(6l - in^v)) sinh {^{6 + in^i^)) ' 

(C25) 



The normalization constant N!^ ^ can be calculated as 



bfc')! bkbi 



the one-particle form factor Fi, is given by equation 
(C14). 

In particular, the pole function for the few lowest 
breathers are 



18 



n: 

sinh(i(6l - mCl) sinh(i(6i + iTrCl) ' 



sinh(i(6l - 2i7rO) sinh(i(6i + 2iO) j 

1 + 2 cos(7r^) 
+ — 

8 cos(7r^) cosh (5 (6* + iirS) ) cosh [ ^ (6* - iirS, 



4 cos 



(^(l+cosK)) 



I r n2 , ^ 

sinh(i(e-^))sinh(i(0 + ^)) V 8cos(^)cosh[i(0 + ^)]cosh[i(0-^)] /' 

^ 1 f 1 

n sinh(i(6l - ifcO) sinh(i(6l + ifcO) V""^ 2 cosh [1(6* + ^tt^) ] cosh [\{0 ~ iir^) ] 



where 



sin(-|-) 
sin(7r^) 



4. Breather-breather form factors of O 

Equation (C9) allows us to express the operator Q in 
terms of exponential operator e"®. Therefore the form 
factors i^b^jj (d) are expressed in terms of the form factors 
of the exponential fields. Taking into account that the 
pole function K^^f^^{9) (C25) goes to zero as a? when 
A: + Z is even and linearly with a when fc + Z is odd, we 
can conclude that only form factors F^^^ (9) with k + I 



being odd are nontrivial. These form factors are written 
as 



where the pole function K^^j^^{9) is given by, 



the normalization factor NP i is determined as 



N, 



bkb, 



= lim A^,f 



a-*0 



bkbi ' 



and the minimal form factor is defined by equation (C7). 
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